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We show that the mean inverse populations of nondecreasing, 
square integrable, continuous-time branching processes decrease to 
zero hke the inverse of their mean population if and only if the initial 
population k is greater than a first threshold mi > 1. If, furthermore, 
k is greater than a second threshold m2 > mi, the normalized mean 
inverse population is at most l/(fc — 1712). We express mi and m2 as 
explicit functionals of the reproducing distribution, we discuss some 
analogues for discrete time branching processes and link these results 
to the behavior of random products involving i.i.d. nonnegative sums. 

Introduction. We consider nondecreasing continuous-time branching pro- 
cesses {zt)t>Q with initial population zq := k> 1, split intensity s{t) and 
offspring distribution p[t) * 5i at time i > 0, where 

p{t) = ^p{i,t)5i. 

i>0 

In other words, zt is the size at time t of a population which evolves as 
follows. During the interval [t,t + dt), each individual living at time t survives 
with probability 1 — s{t) dt + o{dt). For every i > 0, the individual dies and is 
instantaneously replaced by i -|- 1 individuals, independently of the behavior 
of the other individuals living at time t, with probability p{i,t)s{t) dt + 
o{dt). The process {zt)t is pathwise nondecreasing, hence zt > k with full 
probability and the harmonic moments Efc{l/zf} are well defined. 

Harmonic moments of branching processes are needed to build unbiased 
estimators of the offspring distributions from samples of Markov branching 
processes in noncanonical situations, for instance when the state of an indi- 
vidual depends on the number of its siblings. Examples of this situation in 
a discrete-time setting arise in the context of applications to molecular bi- 
ology; see [8] for a presentation and [5, 6, 7] for a mathematical study. Also, 
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harmonic moments reflect mainly lower deviations of the branching process 
and may be viewed as an integrated way to quantify these deviations; see 
[3, 4] and the thorough exposition in [2]. 

In this paper, we provide sharp bounds for Efc{l/zt} if and when E/^.{l/zf} 
decreases to zero roughly like 1/Efc{zt} when t — > oo. The emphasis is on 
explicit and computable upper bounds and the setting is restricted to square 
integrable offspring distributions. 

1. Main results. This section is composed as follows. We first introduce 
some definitions (Section 1.1) and explain a canonical reduction of the model 
(Section 1.2). Then we explain our main results the continuous-time case 
(Section 1.3) and the discrete-time case (Section 1.4). Finally, we sketch 
some relations of our results to previous work (Section 1.5) and describe the 
overall plan of the rest of the paper (Section 1.6). 

1.1. Notation. 

Definition 1.1. For every t > 0, L{t) > denotes a random variable of 
distribution p{t) and 

E{L{t)} = Yip{i,t), M{t):= [\{L{u)]s{u)du. 

^>o 

The function and the quantity h{k) are defined by 

hk{t) := e^'^W Ek{l/zt}, h{k) := sup/ifc(t). 

t>o 

From elementary computations. 

Since the function z 1/z is convex, /ifc(-) is nondecreasing, h{k) is the limit 
of /i/t(-) at infinity, h{k) >l/k = hi^{0) and, for every t > 0, 

e-^^W/fc<IEfc{lM}</i(A;)e-*'^W. 

Hereafter, our aim is to provide sharp upper and lower bounds of h(k). 

1.2. Reduction. In this section, we show that one can assume s(t) = 1 
and L{t) > 1 for every t>0. 

Definition 1.2. For every random variable L such that P(L = 0) / 1, 
L' denotes a random variable distributed like L and conditioned by {L ^0}. 



HARMONIC CONTINUOUS TIME BRANCHING MOMENTS 3 



Definition 1.3. Let (*) denote the condition that the function cj(-) is 
unbounded, where, for every t > 0, 

a{t):= [ {l-p{0,u))s{u)du. 
Jo 

The individuals who survive and those who die and are replaced by ex- 
actly one individual have the same net effect on the overall count zt of the 
population. This proves Proposition 1.4. 

Proposition 1.4 (Canonical reduction). Let (z^)t denote the branching 
process of parameters p'{-) and s'{-) = 1, where p'{t) denotes the distribution 
of L'{a~^{t)). Then the distributions of {zt)t and (^^(t))* coincide. As a 
consequence, provided the condition (*) is met, the quantities h{k) for {zt)t 
and for {z[)t coincide. 

Hence, under condition (*), one can assume without loss of generality (and 
we do assume this from now on) that L{t)>l for every t, almost surely, and 
that s(-) = 1. 

1.3. Results in continuous time. The case when the random tree asso- 
ciated to the branching process is regular is, in a sense, extremal. We first 
deal with this case. 

Theorem 1 (Degenerate case). Assume that p{t) = 5i for every t >0, 
with i>l. Then, for every k <i, h{k) is infinite and for every k>i + l, 

h{k) = l/{k-i). 

In the general case, the parameters mi and m2 defined below play a crucial 
role. 

Definition 1.5. For every integer- valued, not identically zero, square 
integrable L > 0, let 

mi{L} := E{L | L / 0}, m2{L} := K{L^}/E{L}. 

For every L, mi{L} = mi{L'} and m2{L} = m2{L'} (see Definition 1.2), 
and 

1 <mi{L} <m2{L}. 

If L assumes more than one nonzero value, these inequalities are strict. Our 
next result deals with the case when p{t) does not depend on t: 
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Theorem 2 (Homogeneous case). Assume that L[t) is distributed as 
L>0, for every positive t, and that L is square integrable and not identically 
zero. 

(a) For every k < mi{L}, h{k) is infinite. In particular, h{l) is infinite. 

(b) For every k > mi{L}, h{k) is finite and 

h{k) > l/(/e-mi{L}). 

(c) For every k > m.2{Li}, 

l/{k - mi{L}) < h{k) < l/{k - m2{L}). 

Finally, we state our result in the general case: 

Theorem 3 (General case). Assume that there exist finite positive con- 
stants mf and m2 such that, for every t, 

<mi{L{t)} <mt, m2{L(t)} < m^. 

Then the following hold: 

(a) For every k <mi , h{k) is infinite. 

(b) For every k > m^, h{k) is finite. 

(c) For every k > m2 , l/{k - mj") < h{k) < l/{k - 1712"). 

We explain below how to deduce from Theorems 2 and 3 upper bounds 
of h{k) in the ranges mi{L} < k < xn2{L} and < k < m2 , respectively. 
Finally, note that in Theorem 3, mi{L{t)} and m2{L{t)} can be independent 
on t, even when the distribution of L{t) depends on t. 

1.4. Results in discrete time. This section refines results in [7] which we 
recall in Section 2.1, and uses notation from Definition 2.4 in Section 2.1. 

Theorem 4. In the setting of Theorem 2, for every k > m2{L}, 

H{k)<l/{k-m2{L}). 
In the setting of Theorem 3, for every k > m^, 

H{k) < l/{k-mt). 

Overall, the situation is less clear in discrete time than in continuous time; 
note in particular that Theorem 4 does not yield the exact set of populations 
k such that II{k) is finite. 
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1.5. Relation to previous work. In the discrete-time setting, this paper 
deals with offspring distributions di*p such that p{0) 7^ 0, which belong to 
the so-called Schroder's case; see [2], for instance. With the notation of our 
paper, Schroder's exponent a is defined by the relation 

p(0)(l+E{L})" = l. 

When the distribution p is Llog L integrable and when Z{0) = 1, the mar- 
tingale Z{n)/{1 +E{L})"' converges almost surely to a random variable W, 
which is almost surely positive and finite. Introducing i.i.d. copies (VFi)j>i 
of W and using the fact that (1 + E{L})"/Z(n) is a positive submartingale, 
one gets 

Hik)=Eil/{Wi + --- + Wk)). 

Likewise, in our continuous-time setting, when zq = 1, Zt/e^'^^^'> converges 
almost surely to f/;, say, and using i.i.d. copies (it'i)i>i of w, one gets 

h{k)=W.{l/{wi + --- + Wk)). 

Returning to the discrete-time setting, W has distribution dx, the func- 
tion is continuous (see [1]), and there exist positive and finite constants 
ci and C2 such that, when x ^ 0, 

cix^-"<u;(x)<C2X^~°. 

Using this, one can check that Ei(l/iy^) is finite if and only if 

a> [3. 

Likewise, starting from Zq = k,W is distributed like the sum of k i.i.d. copies 
of W starting from one individual. Hence, H{k) is finite if and only if 

ka > 1. 

One sees that, in a loose sense, the initial population k plays the role of 1//?. 
For instance, by convexity, it is a quite general and easily-demonstrated fact 
that if Ei{l/W^/'') is finite, then H{k) = Kk{l/W) is also finite. For more 
remarks on this, see Section 1.5 of [7]. 

1.6. Plan. The rest of the paper is organized as follows. In Section 2, 
we recall some results from [7] about the harmonic moments of branching 
processes in discrete time. We then refine these results and show how to 
apply them to the continuous-time setting. This section culminates with the 
statement of Theorem 5, which we prove later. In Section 3, we prove The- 
orem 1 and state some additional remarks. Section 4 presents a stochastic 
recursion, basic for our study. This allows the exact determination of the 
initial populations k such that h{k) is finite. In Section 5, we prove part 
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(b) of Theorem 2. In Section 6, we prove Theorem 5, stated in Section 2. 
In Section 7, we use the upper bounds for initial populations above m2 to 
deduce upper bounds for smaller initial populations, namely between mi 
and m2. Finally, we explain in Section 8 how our results on branching pro- 
cesses yield the limits of specific random products, which involve some sums 
of nonnegative i.i.d. random variables. 

2. Continuous-time and discrete-time harmonic moments. 

2.1. Discrete-time harmonic moments. We consider branching processes 
in discrete time (Z(n))„>o with reproducing distribution p*Si. Hence, in- 
troducing i.i.d. random variables L{n,i) > 1 with distribution p*6i, we set 

Z[n] 

Z{n+l) := J2 L{ri,i). 

1=1 

As is well known, continuous-time branching processes are limits of discrete- 
time ones. To see this, for every distribution p, we introduce a family {Zu)u 
of branching processes in discrete time, based on p and defined as follows: 

Definition 2.1. For each u in (0, 1], > denotes a random variable 
of distribution pu given by 

(1) Pu := (l - u)5o + up. 

Let {Zu{n))n>o denote a nondecreasing, integrable Galton-Watson pro- 
cess, indexed by the nonnegative integers, starting from Zu{0) := k and with 
reproducing distribution 5i *Pu- Then, when u — > 0, Z„ converges to {zt)t 
in distribution, at least in the sense that, for each fixed t, Zu{\t/u\) con- 
verges to zt in distribution. Thus, results about the harmonic moments of the 
discrete-time branching processes when suitably uniform over the pa- 
rameter u, yield information about the harmonic moments of the continuous- 
time branching process {zt)t- To develop this strategy, we borrow some def- 
initions and results from [7], with some modifications. 

Definition 2.2 ([7]). For every nonnegative m, (70, m 

(s) := s"*. For every 

nonnegative m and every positive c, 

ffc,m(s) :=(l + m(l-s^))-^/^ 

For every nonnegative c, Qc denotes the space of nonnegative random vari- 
ables L such that, for every s in [0, 1], E{s^} < gc,m{s) with m :=E{L}. 
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We recall the following facts and remarks from [7]. First, 

(7c,m(s)=:E{s^-'"} 

is the generating function of a nonnegative random variable Lc^rn^ every 
moment of which is finite, such that E{Lc,m} = and whose distribution is 
either 5m if c = 0, or is supported by the set of nonnegative multiples of c if 
c is positive (viz., the distribution of Lc,m/c is negative binomial). We stress 
that one compares the generating functions of L and of -Lc^m, two random 
variables whose means are equal, but that, the case when c is an integer 
being excluded, Lc,m is not an integer-valued random variable. Hence there 
does not exist, in general, any branching process based on the distribution 
of Lc^m- 

Second, every L in Qc is square integrable and its variance cr^{L} is at most 
c7^{Lc,m} = cm{m + 1). Finally, the family {Qc)c>o is increasing and every 
square integrable L belongs to Qc if c is sufficiently large. Using Definition 2.3 
below, we rephrase this as Proposition 2.6. 

Definition 2.3. For every integer- valued random variable L > 0, let 
e:{L}:=inf{c>0;LGgj. 
When the context is clear, we use £ for '^{L}. 

Our next definition is the analogue in discrete time of Definition 1.1: 

Definition 2.4. For every integer-valued random variable -L > 0, let 
Hn{k) := (1 +E{L})"Efc{l/Z(?i)}, H{k) := Imi^ H^ik) = sup Hn{k). 

Proposition 2.5 ([7]). For every integer k > ^{L}, 

H{k) < l/{k-(t{L}). 

Proposition 2.6 below shows that the result of Proposition 2.5 is not empty 
since ^{L} is finite, at least for (in fact, exactly for) square integrable random 
variables L. 

Proposition 2.6 ([7]). For every non-square-integrable L, is in- 

finite. For every square integrable L, ^{L} is finite and 



E{L}(E{L} + 1)' 
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2.2. From discrete time to continuous time. We apply the results of the 
preceding section to continuous-time branching processes. We recall that 
Definition 2.1 in Section 2.1 introduces random variables such that 
E{Lu} = l + uK{L}. 

Definition 2.7. Let 

Hliik) := (l + uE{L})"Efc{l/Z„(n)}, := lim Hli{k) = supHli{k). 

n^oo „>o 

As a direct application of results in discrete time to our setting, assume 
that pu belongs to Gc{u) ^or every sufficiently small u. Then, for every n, 
every sufficiently small u and every k > c(u). Proposition 2.5 in Section 2.1 
yields 

H^{k) < (1 + uK{L}mi/{Zu{n) - c(n))} < l/{k - c{u)). 
Hence, H^^^^^ {k)<l/{k- c(u)) and 

hk{t) <limJnfl/(/c-c(n)). 
This motivates Definition 2.8 and proves Proposition 2.9 below. 

Definition 2.8. For every integer-valued random variable L>0, let 

c:o{L}:= sup e:{L„}. 

0<«<1 

When the context is clear, we use Co for €o{L}. 

Proposition 2.9. For every integer k > <to{L}, 

h{k) < l/{k - <to{L}). 
Hence, for every nonnegative t, 

Efc{l/zt} < Efc{l/(zj - e:o{^})} < e~^'/{k - U{L}). 

2.3. Identification of Co- The main drawback of the results in Sections 
2.1 and 2.2 is that the parameters ^{L^} (and hence, presumably, the value 
of <tQ{L}) are virtually unknown. In fact, while the exact value of for a 
given random variable L, may indeed be difficult to obtain, it transpires that 
the value of <tQ{L}, which is the only one relevant in the continuous-time 
setting, is quite simple. 

Approaching this result, we first note that is a monotone function 

of the parameter u in (0, 1]. 
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Proposition 2.10. For every c, u and L, if Lu belongs to Qc, then L 
belongs to Qc- In other words, the function u i— > C{Lu} is nonincreasing for 
u in (0,1]. As a consequence, the parameter <t(){L}, defined as a supremum, 
is also given by 

e:o{L} = iime:{L„}. 

Hence ^{L} describes p and ^o{L} describes the limit oi pu when n ^ 0. 
Our next result is anecdotal and stated for the sake of completeness. It 
includes a generalization of Definition 2.4. 

Proposition 2.11. For any positive integer k and any real number b < 
k, let 

hk{t,b) :=e™*Efc{l/(zt -6)}, h{k,b) := sup/ifc(t, 6). 

t>o 

The sequences {h[k))k>i and {h{k,b))k>b CLf^ then convex. 

Our main result in this section provides some tight bounds for ^{L} and 
the exact value of ^q{L}, for every square integrable L. 

Theorem 5. For every integer-valued, nonnegative, square integrable L, 

v^2{L} - < ^L] < m2{L}. 

As a consequence, 

(to{L} = m2{L}. 

The upper and the lower bounds that the first part of Theorem 5 yields 
for (t{Lu} both converge to m2{L} when 0, hence the second part is a 
direct consequence of the first part. Finally, Theorem 5 and Proposition 2.9 
imply part (c) of Theorems 2 and 3. Together with Proposition 2.10, they 
also yield Theorem 4. We prove Theorem 5 in Section 6. 

2.4. Proofs. 

Proof of Proposition 2.10. Since E{s-^"} = i - n + uE{s^}, it is 
enough to show that gc,um. <1 — u + ugc,m.- In turn, this follows from the 
convexity of the function m gc,m{s) for every fixed s. □ 

Proof of Proposition 2.11. This is elementary. For every positive 
i.i.d. Q and Ci and every k>l, let (T^ := Ci + • • • + C,^ and rjk := E{l/(Tfc}. 
Then the sequence {r]k)k is convex. To see this, note that the convexity 
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is equivalent to {rjk — rjk+ijk being nonincreasing and that this holds true 
because 

% -%+i =IE{C/((C + crfc)o"fe)} 

and because {ak)k is nondecreasing. Applying this convexity to Zu{n) and 
considering the limit when n — > oo shows that is a convex sequence for 
every u. Taking the limit when n ^ shows that the sequence h is convex. 
The proof for h) is similar and hence omitted. □ 

3. Regular trees. The simplest case is when the trees are regular, that is, 
when the branching process is ruled by a distribution p : = 5i with i>l. We 
already mentioned that the case (1 — X)6q + X5i with A in (0, 1] is equivalent, 
since this corresponds to a change of time of the branching mechanism. The 
associated random tree is regular with degree i + 1 and the length of every 
edge is exponential of mean 1/A. The distribution pu is Bernoulli or z with 
Puii) =uX. Such Bernoulli distributions belong to Qi. In fact (see [7]), 

^0 = mi = m2 = i. 

Hence, for every k>i + l, 

h{k) < l/{k-i) 

and this method cannot yield a finite upper bound in the k = i case. This is 
for good reason since, as we show below, h{i) is in fact infinite. 

The distribution of zt when zq = 1 is an affine transform of the negative 
binomial distribution of index i and mean e'*'*. More precisely, for every 
\v\ < 1, 



Hence, for every positive t, 

z(e^' — 1) 

and hi increases from hi{0) = 1/i to h(i) = oo. This shows that is un- 
bounded when k < i. On the other hand, for every k > i + 1, the result 
above is sharp since (see below) the function hk increases from hk{0) = 1/k 
to 

h{k) = l/{k-i). 

For branching processes in discrete time (see [7]), the true value of ^{L} is 

^^"^^ log(l - A) ' 
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while Proposition 2.6 yields the a priori bound 

The function A i-^ decreases from (t = (to = i~ at A = O"*" to (t = 0^ at 
A = 1~. Hence, this function is discontinuous at A = since A = yields 
linear trees such that (L = 0. 



4. A stochastic recursion. To deal with the finiteness of h in the general 
continuous case, we return to the recursion in discrete time. Recall that L' 
denotes a random variable distributed like L, conditioned by {L 7^ 0}, and 
that 

mi{L} := E{L | L / 0}, m2{L} := E{L'^}/E{L}. 
Hence, for every u in (0, 1], 

mi{L} = mi{L'} = mi{L„}, m2{L} = m2{L'} = m2{L„}. 

4.1. Results. Our first aim in this section is to prove Proposition 4.1. 

Proposition 4.1. Assume that L is such that ^o{L} is finite. Then 
h{k) is finite if and only if k > mi{L}. 

Corollary 4.2 follows directly and is anecdotal. 

Corollary 4.2. For any L such that Co is finite, h{k,b) is finite if and 
only if k > max{mi, b}. 

Our second aim is to characterize the sequence {h{k))k as solution of a 
stochastic recursion. 

Proposition 4.3. For every L such that G!o{-^} finite and every k > 
mi{L}, 

(2) {k-E{L})h{k)=kE{h{k + L)}. 

This is equivalent to the fact that, for every k > mi{L}, 

(A; - mi{L])h{k) = kE{h{k + L')} = kE{h{k + L)\L^ 0}. 
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4.2. Remark. As was expected, the distributions p and pu yield the same 
functional h. To check this a posteriori, assume that h solves the recursion 
in Proposition 4.3 for L of distribution p and that /i" solves the recursion 
in Proposition 4.3 for (Caution: corresponds to the continuous-time 
process with offspring distribution pu and to the discrete-time process 
Zu-) Then = uK{L} = um and Proposition 4.3 for becomes 

{k - um)h}'{k) = k¥.{h}'{k + Lu)}. 

By definition of L^, 

E{/i"(A; + Lu)} = (1 - u)h''{k) + vE{W{k + L)}. 

Rearranging and dividing both sides by u, one gets 

{k - m)h}'{k) = kEih'^ik + L)}. 

Since h and /i" satisfy the same relation, h = h^ for every u in (0, 1]. Finally, 
since the distributions of and (i-p(o)) coincide, L and L' also yield 
the same function h. 



4.3. Proofs. 

Proof of Proposition 4.1. Assume without loss of generality that 
L > 1 almost surely and let m := mi{L} = E{L}. Every H^{k) is finite, 
-ffg (fc) = l/k, every sequence {H^{k))n is nondecreasing and, for every n > 1, 

(3) H^iik) = (1 + 7xm)Efc{//^i(Z4l))}. 

First, since Zu{l) = k with probability (1 - u)^ on {Z„(0) = k}, 

K{k)>{l + um){l-ufHl^^{k). 

If k < m, then (1 -|- um) (1 — u)^ > 1 for every u sufficiently small. This 
implies that H'^{k) is infinite for every u sufficiently small, hence h{k) is 
infinite. 

Herefar, we assume that k > m. Since is nonincr easing, one has 

H^_-^{Zu{l)) < Hii_^{k + 1) on the event {Z„(l) >k + l}. The probability 
of {Zu{l) > A; -I- 1} when Z„(0) = A; is 1 - (1 - u)^, hence (3) implies that 

H:{k) < (1 + nm)((l - u)'H:^,{k) + (1 - (1 - u)')H:^,{k + 1)). 

Since 1 — (1 — u)^ < ku and < -ff", we have 

K{k) < (1 + 'um)((l - u)''H^{k) + + 1)). 

Since k> m, for every positive and sufficiently small u, (1 -|- um){l — u)^ < 1. 
Hence, for u sufficiently small, 

ku{l + um) 



H:i{k)<H:{k+i) 



1 — (1 + um){l — u)^ 
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Letting n oo, the same inequality holds true for H^{k) and H'^{k + 1). 
When n — > 0, for every fixed k > m, the fraction on the right-hand side 
converges to the finite limit k/{k — m), hence if h{k + 1) is finite, then h{k) 
is finite too, and 

h{k)<h{k + l)k/{k-m). 

Finally, if h{k) is finite for k sufficiently large, and indeed h{k) is finite for 
k > (tQ{L}, then h{k) is finite for every k> m. 

The case A; = m is similar to the case k <m. For every k, 

K{k) > (1 + Mm)((l - u)^Hl{k) + ku{l - u)^-^Hl{k + 1)). 

One concludes by writing this as H^{k) > a'^{k)H^(k + 1) for an explicit 
o5^(fc) and by checking that a'^{m) oo when ti — > 0. □ 

Proof of Proposition 4.3. The sequence {H'^{k))k solves the recur- 
sion 

H^{k) = {l + um)¥.k{H^{Zu{l))}. 

When Zu{0) = k and n — > 0, either (i) every individual survives, this hap- 
pens with probability 1 — ku + o{u) and then Zu{l) = k, or (ii) exactly one 
individual dies and is replaced by 1-l-L individuals, this happens with prob- 
ability ku + o{u) and then Zu{l) = k + L. The probability that at least two 
individuals die is o{u). Hence, 

H''{k) = (1 + um){{l - ku)H''{k) + kuW^iH^'ik + L)}) + o(n). 

It follows that h{k) = \\m.u^{) {k) , if finite, solves the recursion stated 
in the proposition. The finiteness of h{k) is crucial in this argument since 
one must cancel h{k) on both sides of the equation, and this finiteness is 
guaranteed by Proposition 4.1. □ 

In relation to the last sentence of the proof of Proposition 4.3, note that 
if IE{L} < < mi, equation (2) in Proposition 4.3 is still valid in the sense 
that both sides are infinite. 



5. The first threshold. Proposition 4.3 in Section 4.1 shows that h{h) 
cannot be finite if fc < mi. Hereafter, we assume that L >1 almost surely 
and that k > mi = m := E{L}. In the special case p = 6i, L = i almost surely 
and iterating Proposition 4.3 yields 

{k - ^)h{k) = {k+ ji)h{k + {j + > ^ . 

for every integer j > 0, hence h{k) > l/{k — i). Since the reversed inequality 
was established above, the proof of Theorem 1 is complete. 
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Returning to the general case, we now prove part (b) of Theorem 2, which 
can be viewed as an extremahty property of the deterministic branching 
mechanisms since it asserts that, for every other branching mechanism, for 
every k > mi{L}, 

h{k) > l/{k-mi{L}). 

Proof of part (b) of Theorem 2. One can assume without loss 
of generality that L > 1 almost surely. Let m := E{L}. Since h{k + L) > 
l/{k + L), we have 

{k - m)h{k) > kE{l/{k + L)} > k/{k + m). 

Thus, the a priori lower bound h{k) >l/k can be improved to the lower 
bound 

k 



h{k) > 



{k — m){k + m) 



One can iterate this reasoning as follows. Assume that h[k) > R{k), where 
R{y) denotes a rational fraction, convex on y > m. Then 

E{h{k + L)} > E{R{k + L)} > R{k + m), 

and hence h{k) > T{R){k) for a new rational fraction T{R), convex on y > m 
and defined as 

y — m 

Starting from R{y) := 1/y and iterating this, one gets h{k) > T'^{R){k) with 



n+l\y + nm y — m, 

Letting n — > oo, this yields h{k) > l/{k — m). Finally, h{k) > l/{k — m) 
as soon as some strict convexity is involved, that is, for every distribution 
P 7^ 5m- For instance, one can write [k — m)h{k) as 

kE{h{k + L)} > kE{l/{k + L - m)} > k/{k + E{L} - m) = 1. □ 

6. The second threshold. Our aim in this section is to prove Theorem 5 
stated in Section 2.3. 

6.1. Results. It seems that the parameter m2 introduced in Definition 1.5 
enters the picture mainly because of the technical result in Lemma 6.2 below. 



Definition 6.1. For every positive c and every s in [0, 1], let 
(4) ifcis) := c(l - E{s^}) - E{L}(1 - s''). 
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Lemma 6.2. Let c> 0. Then c > m2 if and only if fds) > for every s 
in the interval [0, 1] . 

Lemma 6.3. For any given L, ^o{L] >m2{L}. 

Lemma 6.4. For every L, <t{L} <m2{L}. Since m2{L'} = m2{L} and, 
for every u in (0, 1], m2{Lu} = m2{L}, this implies that <t{L'} < m2{L} and 
e:{L„}<m2{L}. 

This implies the upper bounds in Theorem 5. 
Finally, we note that for every L and every u in (0, 1] , 

iL{L'} < iL{L} < < ito{L} = m2{M}, 

where M may be any of the random variables L' , L or 

6.2. Proofs. 

Proof of Lemma 6.2. First, ipc{s) ='i/'c(s'^) with 

Ms) ■■= c(l - E{s^/^}) - E{L}(1 - s). 

Then i^d^) = V'c(l) = ^i^d "0" (1) = ]E{L}(1 — m2/c). Hence, for any positive 
c < m2, some values of (fds) are negative. 

Conversely, assume that c > ro2. If V^c is convex, then (/^c > on [0,1]. 
To show that ^l^c is indeed convex, we first note that ^l^'^{s^^^) has the same 
sign as E{L(c — L)s^}, which increases as c increases. Thus, it is enough to 
show that this last expectation is nonnegative when c = m2. In turn, this 
condition is equivalent to X ^ on [0, 1], where 

Xis) := E{L'^}E{Ls^} - E{L}E{L'^s^}. 

One sees that xi^) = E{K{s)}/2, where 

i^(s):=LiL2(Li-L2)(s^^-s^^) 

and Li and L2 are two independent copies of L. Since the function s^ 
is nonincreasing, K{s) > almost surely and hence x{s) > 0. This concludes 
the proof of Lemma 6.2. □ 

First proof of Lemma 6.3. For every positive c> G^o, Lu belongs to 
Qc for every u, in particular when u — > [we recall that the family {Qdc 
is increasing]. Rearranging both sides of the inequality which characterizes 
the property that belongs to Qc in Definition 2.2 and taking logarithms, 
one sees that this inequality is equivalent to the fact that u;l(s,?x,c) < for 
every s in [0, 1] where 

(5) wl(s, u, c) := log(l + uE{L}{l - s^)) + clog(l - n(l - E{s^})). 
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Since ujl{s,u,c) = —ipc{s)u + o{u) when u ^ 0, one gets ^c{s) > for every 
s in [0, 1], that is, c > m2; see Lemma 6.2. Finahy, every c > Co is such that 
c > m2, hence Lemma 6.3 holds. □ 

Second proof of Lemma 6.3. Theorem C in [7] states that for any 
L in Qc, 



We apply this to L„. Since E{L„} = uK{L} and K{Ll} = uE{L'^}, this 
inequality is equivalent to €{Lu} > m2{L}(l + o(l)) when u ^ 0, hence 
<Lo{L}>m2{L}. □ 

Proof of Lemma 6.4. Using the notation from the proof of Lemma 6.3, 
we must show that u}l{s,u, xn2) < 0. 

Since log(l + u) <v for every v, u!l{s,u,c) <<^l{s,u,c) for any c with 



The first term in ?l(s, -,0) is linear with respect to u. The second term has 
the form /31og(l — au) with positive a and /3, hence it is concave with respect 
to u. Thus, <^l(s,-,c) is concave. Since ?l(s,0,c) = 0, we have (^l{s,u,c) <0 
for every u in [0,1] if and only if the first derivative of (;l{s,-,c) at is 
nonpositive. Since this derivative is —ipc{s), Lemma 6.2 above shows that 
c) < if and only if c > rri2. □ 

7. Between the two thresholds. Theorem 2 describes h{k) when k <mi 
and when A; > m2. In this section, we study the regime mi < A; < m2. We start 
from the observation, drawn from Proposition 4.3, that for any A; > mi. 



One can iterate this trick to deal with values of A; > mi such that A; + 2 > m2 , 
and so on. This yields Proposition 7.1. 

Proposition 7.1. Let A; > mi and let i denote any nonnegative integer 
such that k + £> m2 . Then 



a^{L}<cK{L}{E{L} + l). 



u, c) := uE{L}{l - s^) + clog(l - u(l - E{s^})). 



{k-mi)h{k) <kh{k + l). 
For instance, if m2 — 1 < A; < m2, then A; + 1 > m2 and 

h{k + l) < l/(A: + l-m2). 
Thus, for every A; > mi such that A; + 1 > m2, 

h{k) < k/{{k - mi)(A; + 1 - m2)). 
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Every denominator in the upper bound for h{k) is positive. A strange 
feature of this result is that the best upper bound for h{k) does not al- 
ways correspond to the smallest value of i such that k + i > m2- Simple 
computations show that the optimal choice is the smallest value of £ such 
that 

k + l> mi = m2 H — > m2. 

mi — 1 mi — 1 

8. Branching processes and i.i.d. sums. 

8.1. Random products. Our starting point in this section is the recursion 
of Proposition 4.3 which motivates the following definition: 

Definition 8.1. Let X and Xi denote i.i.d. nonnegative integrable ran- 
dom variables. Let 5*0 := and Sn '■= Xi^ \-Xn for every n>l. For every 

positive X, let 

Rn{x,X) :=[j(l + M^l g^j^^x) ■.= limsupE{Rnix,X)/n}. 

We first deal with the finiteness of the functional q, then we relate it to 
the harmonic setting. Lemma 8.2 lists some elementary facts, the proofs of 
which are omitted. 



Lemma 8.2. Let X denote a nonnegative, nonzero random variable. 

(i) The function x ^ q{x,X) is nonincreasing. 

(ii) For every nonnegative n, ¥.{Rn{x,X)} > 1 + (n -|- l)E{X}/x. 

(iii) Hence, g{x,X) >E,{X}/x. In particular, q{x,X) is always positive. 

(iv) For every positive X, g{Xx,XX) = q{x,X). 

8.2. Results. 

Definition 8.3. For every nonnegative, not identically zero, integrable 
X, let 

mo{X} := F{X = 0}mi{X} = mi{X} - E{X}. 

Note that mo{X} satisfies 

(1 +E{X}/mo{X})P{X = 0} = 1. 

Proposition 8.4. For every nonnegative integrable X, q{x,X) is infi- 
nite for every x < mo{X}. 
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Proposition 8.5 implies a dichotomy result, stated as Corollary 8.6. 

Proposition 8.5. There exist finite positive functions v{-,X) andw{-,X) 
such that, for every positive x > mo{X} and every positive y, 

q{x + y,X)< g{x, X) < q{x + y, X)v{x, X)w{x, Xf . 

Corollary 8.6. Either q{x,X) is infinite for every x > mo{^}, or 
g{x,X) is positive and finite for every x > mo{X}. 

Theorem 6 links some functions g to the harmonic moments of branching 
processes: 

Theorem 6. For every square integrable L with nonnegative integer 
values and for every integer k > mi{L}, the normalized limiting harmonic 
moment h{k) which corresponds to the branching mechanism based on L 
satisfies the relation 

_ g{k-E{L},L) _ g{k-mi{L},L') 
E{L} ~ mi{L} 

Corollary 8.7 below is a direct consequence and complements Corollary 8.6. 

Corollary 8.7. LetX denote a nonnegative, integer, square integrable 
random variable. Then g{x,X) is finite if and only if x > mo{-^}. Fur- 
thermore, xg{x,X) — > E{X} when x — > oo. More precisely, for every integer 
x>mo{X}, 

E{X}/x < g{x - E{X},X) < E{X}/{x - m2{X}). 

From assertion (iv) in Lemma 8.2, the same conclusion holds if X is a 
multiple of a nonnegative, integer, square integrable random variable. The 
hypothesis that X is lattice and that x belongs to this lattice might be 
unnecessary. 

Note that one uses the square integrability of X only to ensure the con- 
dition that limsup^j^oo < 1 [which is equivalent to limk^oo kh{k) = 1]. 

For every positive x and every nonnegative random variable X , the func- 
tion ip := g{-,X) solves the recursion 

(6) xip{x)=E{ip{x + X)}{x + E{X}). 

Recall that g{x,X) > E{X}/x in full generality. A consequence of (6) is that 

xg{x,X) >E{X}E{x + X}E{l/{x + X)}, 

hence g{x,X) > E{X}/x in full generality. 
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8.3. Proofs. 

Proof of Proposition 8.4. Let r :=¥{X = 0}. Since Xi is with 
probability r or distributed like X' with probability 1 — r, we have 

E{Rn+iix,X)} 

(7) 

= {l + E{X}/x){rE{Rn{x,X)} + {l-r)E{Rn{x + X',X)}). 

The second part of Lemma 8.2 yields 

E{Rnix + X',X)}>l + {n + l)E{X}E{l/(x + X')}. 

If (1 + E{X}/x)r > 1, this implies that there exists a positive constant c 
such that 

E{Rn+i{x, X)} > E{Rn{x, X)} + c{n + 1). 

Hence, E{Rnix,X)} > and g{x,X) is infinite. Finally, the condition 

(1 + E{X}/x)r > 1 can be rewritten as 

x<E{X}r/{l-r)=mo{X}. 

This proves Proposition 8.4. □ 

Proof of Proposition 8.5. Let T denote the hitting time of the level 
y by the nondecreasing process (5'j)j>o, that is, 

T = mf{i>0;Si>y}. 

Since 5j <?/ if i<T— 1 and Si > y if i > T , 

E{Rn (x, X)\T} < RT+n{x, X) < Rt-i{x, X)R'^{x + St, X) 

and hence 

E{Rn{x,X)\T} < {l + E{X}/x)^R'^{x + y,X), 

where denotes a copy of independent of T. Hence, 

q{x,X) <Q{x + y,X)E{s^}, s:=l+E{X}/x. 

To bound the last term, we use the fact that {T > n} = {Sn < y} and hence 
an exponential inequality yields 

P{T > n} < ii~^E{u^}" 

for every positive u < 1. When ti — > 0, E{n"'*"} —>■ ¥{X = 0}, hence the as- 
sumption on X implies that there exists u such that sEjii'''-} < 1. Hence, 

E{s^} < >n-l}<u~^Yl isE{u^}f 

n>l n>0 
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and the last series converges. This proves the proposition. □ 

Proof of Theorem 6. The choice X = L' gives E{X} = mi{L}. Hence, 
Proposition 4.3 yields that, for every k > mi{L} and every n > 0, 

h{k) = E{Rnix, L')h{k + 5„+i)}, x:=k- mi{L}. 

Since Rn{x,L') and h{k + Sn+i) are nonincreasing functionals of the se- 
quence {Xi)i, a coupling inequality reads 

h{k) > E{Rn{x, L')/n}nE{h{k + 5„+i)}. 

Now, by convexity, 

E{h{k + Sn+i)} > E{l/{k + Sn+i)} > l/{k + mi{L}(n + 1)). 

When n oo, we obtain h{k) > q{x, L')/mi{L}. 

On the other hand, for every nonnegative integer i, 

h{k) < E{Rn{x, L')}h{k + E{Rn{x, L')- Sn < i}. 

First, when L' is square integrable, h{k + i) ~ 1/i when i oo. Second, we 
use the following special case of Cramer large deviations bounds: 

Lemma 8.8. Let X and Xi denote i.i.d. nonnegative random variables. 
For every y < E{X}, there exists r < 1 such that for every n, 

¥{Xi + • • • + X„ < ny} < r". 

We apply this lemma to X = L' . This yields P{5„ < ny} < r" with r < 1 
for every y < mi{L}. Third, since L' > 1, Rn{x,L') is at most the product 
of 1 + mi{L} / (x + j) from j = to n, which is 0{n^^^^^). Thus, for every 
y < mi{L}, using the above for i := [ny\ , one gets 

h{k) < Q{x,L'){n/[ny\+o{l))+0{n"''^^^)r''. 

This suffices to prove Theorem 6. □ 

First proof of Lemma 8.8. Let Xj := mi{Xi,t}. For any t, 

Xi + ---+Xi<Xi + --- + Xn 

and, if t is sufficiently large, y < E{X|}. Since X* is almost surely bounded 
and the sequence {X^)n>i is i.i.d., standard large deviations bounds apply 
and hence there exists r < 1 with 

P{Xi + --- + Xn<ny}< F{Xi + • • • + X* < ny} < r". □ 

A second proof of Lemma 8.8 uses the fact that for any nonnegative 
integrable random variable X, E{e~^^} = l — tE{X} + o{t) when t is positive 
and t — > 0. 
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Second proof of Lemma 8.8. Assume first that X is integrable. By 
Chebyshev's inequality, the result of the lemma holds with 

r:=e*^'E{e-*^} 

for every positive t. Since X is integrable and positive, its Laplace transform 
E{e-*-^} is differentiable at t = 0+, that is, E{e-*-^} = 1 - fE{X} + o{t) 
when t is positive and t Since e*^ = 1 + + o{t) and y < K{X}, one 
obtains a value r < 1 for any t positive and sufficiently small. 

If X is not integrable, truncating X and every X„ as in our first proof 
yields the result. □ 
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